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Abstract 

Critical dynamics of the Nambu- Jona-Lasinio model, coupled to a constant 
electromagnetic field in D = 2, 3, and 4, is reconsidered from a viewpoint 
of infrared behavior and vacuum instability. The latter is associated with 
constant electric fields and cannot be avoidable in the nonperturbative frame- 
work obtained through the proper time method. As for magnetic fields, an 
infrared cut-off is essential to investigate the critical phenomena. The result 
reconfirms the fact that the critical coupling in D = 3 and 4 goes to zero even 
under an infinitesimal magnetic field. There also shows that a non-vanishing 
F^uF'^'^ causes instability. A perturbation with respect to external fields is 
adopted to investigate critical quantities, but the resultant asymptotic expan- 
sion excellently matches with the exact value. 
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I. INTRODUCTION 



The four-fermi interaction model by Nambu and Jona-Lasinio (NJL) has been dis- 
cussed to investigate the dynamical symmetry breaking (DSB) in a number of cases in two, 
three, and four dimensions. Especially interesting situations are found such that NJL is 
coupled to external sources, which enables us to peep into detailed structures of DSB and 
obtain information on the chiral symmetry breaking (xSB) in the QCD vacuum, the planar 
(2 + 1-dimensional) dynamics in solid state physics, or the early universe when coupled to a 
curved space-time 0. In this respect the NJL model minimally coupled to the electromag- 
netic fields is discussed by many authors to explain how ^SB could be changed under the 
influence of the electromagnetic fields: Klevansky and Lemmer [Q] find that a pure electric 
field opposes xSB to restore chiral symmetry, meanwhile a pure magnetic field enhances 
xSB. The former result has been generalized to non-abelian gauge fields by Suganuma and 
Tatsumi, where they argue about chiral symmetry restoration by (color) electric fields [Q. 
Meanwhile the latter result is further investigated by Gusynin, Miransky, and Shovkovy 
who find that there occurs the mass generation even at the weakest attractive interaction in 
2 + 1-dimension (xSB in 2 + 1-dimension can be realized as a flavor symmetry breaking 
by introducing an additional fermion 0) and in 3 -|- 1-dimension and emphasize it by 
means of the dimensional reduction. This implies that the critical coupling is zero even if 
the applying magnetic field is infinitesimal, which might however contradict with a naive 
consideration. The motivation for this work lies here. 

Our strategy to this issue is: to start with the Euclidean path integral representation 
of the NJL model minimally coupled to constant electromagnetic fields. Following the 
standard procedure, that is, introducing auxiliary fields and integrating with respect to the 
fermion field, we arrive at the pure bosonic path integral consisting of the auxiliary fields 
and a functional determinant. We then rely on the semi-classical approximation and employ 
the Fock-Schwinger proper time method [§| in order to define a functional determinant. 
Although the proper time method can automatically provide a gauge-invariant ultraviolet 
(UV) regularization in terms of the gamma or zeta function similar to the dimensional 
regularization, we introduce a cut-off to grasp physical situations better, which, however, 
still be gauge invariant as well as Lorentz covariant. Moreover an infrared (IR) cut-off e 
must also be introduced, since infrared divergences arise when external fields are coupled to 
a massless state. So far a very little care has been devoted to this fact but it is inevitable in 
order to discuss effects of electromagnetic fields to the NJL model; since we are interested 
in the transition from massless to massive states or vice versa under an influence of external 
fields. From these procedures we obtain the effective potential. 

Another issue comes up when an electric field comes into play. The effective potential 
becomes complex, whose imaginary part implies the vacuum instability because of the cre- 
ation of fermion and antifermion pairs. (The phenomenon is closely related to the Klein 
paradox in the one body Dirac equation [^].) However if we notice that the imaginary part 
behaves as e~"^ with E being the magnitude of electric field, it is negligible when 

E<m^, (1.1) 

which leads us to the situation that the magnitude of electric fields must always be infinites- 
imal. However it is enough to live in this world, since we are interested in a change of a 
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vacuum given by the NJL model to that with external fields, so that we can consider them 
so small as perturbations. With these spirits, we find that there is no notion of criticality 
in the pure electric field case, since that is defined through a transition from a massive to a 
massless state but there is not any state with a mass less than a magnitude of the electric 



field owing to (|1 . 1| ) . On the contrary, in a pure magnetic field case, we find critical couplings. 
They are nonvanishing as far as the infrared cut-off e is kept finite but eventually become 
zero when e — > 0. This reconfirms the results of Gusynin, Miransky, and Shovkovy who have 
relied on the argument in terms of the dimensional reduction. These results are obtained 
through perturbations with respect to external fields but a resultant asymptotic expansion 
is found well matched with the exact value even if the first few terms are adopted. This 
observation would be interesting. 

The paper is organized as follows: in Sec. we develop a general formalism for com- 
puting an effective potential under constant electromagnetic fields. In the subsequent Sec. 
pn] , [TV] , and |V|, a pure electric, magnetic, and a non- vanishing F^^F'^'' (= e^''^" F^yFp„/2) in 
D = A cases are discussed in order. The final Sec. is devoted to conclusion. In Appendix 
^ calculations of the trace and the determinant in deriving the effective potential are given. 



II. GENERAL FORMULATION 

In this section we describe the model and develop a general formulation for obtaining 
the one-loop effective potential of the D-dimensional NJL model. The Lagrangian for the 
NJL model minimally coupled to external electromagnetic fields is 



D = 2A 
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D = 3 



where the electromagnetic coupling constant has been absorbed into the definition of A^. 
Apart from a usual 4-dimensional case, ip is a two-component spinor with gamma matrices 

7m = 0-/. , 75 = -«7i72 , CT/.!. = [7m) 7i/]/2i = e^^ag ; /i = 1, 2 , (2.2) 

in 2-dimension. For the 3-dimensional case, a spinorial representation of the Lorentz group 
is given by two-component spinors, so that corresponding gamma matrices are 2x2. There 
is no chiral symmetry. In order to be able to discuss chiral symmetry, we introduce another 
fiavor such that |^ 

^=(^j), ^ ^ ^2) ^ {4^3 4^3) , (2.3) 
and 4x4 gamma matrices 

7.= (^0^ _°J ;/i = l~3, 74=(;j), 75 = 71727374 =( ) . (2.4) 

The chiral symmetry realizes as 

^ — , e*""^*^ , ^ — ^ e*^^5^ , (2.5) 
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yielding a global U{2) symmetry which is broken by a mass term into U{1) x U{1). 
The partition function of the model is read as 



Z[A] = JVipViljexp 
IVaVTzVipVip exp 



d^xC 



d^x 



29' 



(2.6) 



where auxiliary fields, a and tt, have been introduced as usual; 

^.P=|^75 for D = 2,4, 

\ vri74 + 7r275 for D = 3 . 

The fermionic integration gives the functional determinant: 



(2.7) 



Det 



^^{d^-iA^) + {a + iTz ■ F) 



(2.8) 



We then perform the semi classical approximation, that is, shift, cr — m + a', tt ^ tt', and 
assign a' and tt' as the new integration variables to find 



Z[A] = exp [-VTVoim)] (l + 0(2-loop) 



with 



Vn(m) = 



m 
2? 



VT 



In Det 



lfi{d^ - iAf,) + m 



(2.9) 



(2.10) 



where V is the (D — l)-dimensional volume of the system and T is the Euclidean time interval. 
(This semiclassical approximation would be more justified by introducing fermion pieces 
and taking N ^ 00.) It should be understood that the terms of 0(2-loop) are given by 
integrations with respect to a' and tt'. Vo{fn) in ( p.lO| ) is the one-loop effective potential 
from which we can see the phase structure of the model. 

In order to make (12. 81) well defined, we normalize it as follows: 



Det [7^ ((9^ - iA^) +m] 



Det [7^9^] 

IYLn(7^(9^ - tA^) + m) - TrLn(7^9, 
iTrLn(^-(9^ - iA^,f - ]-a^yF^^ + 



-TrLn(-9; 



(2.11) 



where the trace operation, designated by Tr, must be taken with respect to the space-time 
as well as the gamma matrices, whereas tr implies that only for the gamma matrices. With 
the use of the identity. 



In if = — lim 



drr'-^e-^^ 



-Vis] 



(2.12) 
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If) can be rewritten as 

2s-oio ^ 



1 /■°° 

- lim / drr''^^ 

2 s^o Jo 



-trl jd^x{x\e^^^x) 



(2.13) 



where, in the second hne, F^^ has been assumed to be constant and a^^ = [^^,'~^y\/2i. The 
kernel (xje'^"'^^"*^''-' ]|a;') with a constant F^^ can be calculated by means of the proper time 
method [Hi as 



(X e 



x') 



exp 



(47rr)^/2 
X exp 



det 



/ sin(rF) ^ 



- — {x - x')^{tF cot tF)^^{x - x') 
4r 



(2.14) 



where F denotes a. D x D matrix whose components are F^^. Combining (|2.13|) and (p.l4|) 
we obtain 



trl 



^-0 2(47r)^/2 Jo 



/■oo n r 9 

dTT'---^ {e^"™ GoirF) - ij 



(2.15) 



where GoirF) reads 
Gd{tF) 



tr (e'S'"^'"'"'^'"^ ^ 
trl 



det 



/ sin(rF) 



tFd coth(rF^) for D = 2,3 

T^F+F_ coth(rF+) coth(rF_) for D = A , 



(2.16) 



with 



F2 

F3 



E 



(2.17) 



F± = {\B + E\±\B - E\] /2 . 



(Details are shown in Appendix 0.) As was mentioned above trl in ( p.l5| ) and ( |2.16| ) is 
2DI2 D = 2,A and 2(^+^)/2 = 3. Although the integral ( ^TTsD has entirely been 

regularized if an analytic continuation is made for s, in order to grasp a physical situation 
better, an UV cut-off A^ is introduced as is done in the ordinary gap equation |l| . Moreover, 
an IR cut-off e is indispensable, since if m? = and F ^ 0, the integral ( |2.15| ) becomes 



divergent. Therefore we consider instead of ( p.l5|) 

trl 



2(47r)^/2 Ji/A 



r dTT-^-^ fe-"(™'+^)GB(rF) - l} . 

Jl/A2 ^ ^ ^ J 



(2.18) 



5 



This regularized Id^ enables us to investigate dynamics near the massless region, so does 
that of xSB. Now the well-defined effective potential is given by 

^2 trl 



m 



2^2 2(47r)^/2 yi/A2 



j^^dTT-^-^ {e-^('™'+^)Gz)(rF) - 1} . 



(2.19) 



We then introduce dimensionless quantities, x = (m^ + e)/ Is} and T = F j {^nu = 
F^jy/A^), obeying e/A^ < a; < 1 and \ J-'fiu\ < 1- (Recall that mass dimension of the gauge field 
is always one because of inclusion of the coupling constant.) Therefore the (dimensionless) 
effective potential Vd{x) is read as 



2(47r 



\D/2 



trl A^ 



(47r 



.D/2 



trl g^A^~^ 



X 



A2 



(2.20) 



The stationary condition, dV ^{m) / dm = 0, is 

(47r)^/2 



dVDix) _ 2m 

X 



d{m/A) 
where foix) is defined as 

foix) 



trl g^A^-^ 



+ fnix) 



d roc 

dx 



/CO n 
rfrr---^e-™Gz)(r^) 



(2.21) 



(2.22) 



Thus the equation for non-trivial solutions of ( p.21|) , the gap equation, is 

f47r)^/2 

Hereafter we designate the stationary point, the solution of the gap equation, as x*(= 
(m*2 -I- e)/A2). The stability condition d"^!) oijn) / dm^\m* > 0, gives 



d^Vnix 



d{m/Ay 



4m" 
X" 



foix*) > 



(2.24) 



where use has been made of ( |2.23| ), and then, the absolute minimum condition, Voifn*) 
VDiO) < 0, leads to 



Voix*) - VD{t/A^) 



D/2 



/A2 ^ ' trl^2^^^2 ^2 



e/A- 



dx [/d(x)-/d(x*)]<0. 



Finally, we list the integral expression for the effective potential in each dimension. 



V2{X) 

Vsix) 
V,{x) 



^{x- — ]+ I dTr-^e-^%rJ^2) coth(r^2) , 



2n 

9^ 
27r3/2 

C/2A 

47r2 



e 

A2 



X - — ) + / c/rr-2e-"^(rj^3) coth(rJ3) 



X 



A2 



/oo 
dr r-=^e-™(r2j-+J^^„) coth(rJ^+) coth(rJ^_ 



(2.25) 

(2.26) 
(2.27) 
(2.28) 



where we have ignored the terms independent of mass and fields. 
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III. CONSTANT ELECTRIC FIELD 



In this section analyses of a constant electric field are made. (The case reads, covari- 
antly in the Minkowski metric, F^^F^"' < as well as Ff,^F^"' = e^^xpF''" F^p /2 = in 
4-dimension.) As was mentioned in the introduction, an imaginary part arises in the effec- 
tive potential: the integral in ( p.22|) becomes 



dTT~--^e-^'-'{rS)coth{T£) for D = 2,3,4, (3.1) 

where S is the magnitude of the electric field in the A^ unit (0 < £^ < 1). To see the 
imaginary part, go back to the Minkowski space via an analytic continuation, S —iS, 
obtaining 

/oo n 
dr T-T-^e~^%r£) cothirS) 

—^PJ^ drr---h~^''{TS)cot{TS)-i7rY,^i^-^j e'^r- (3.2) 

with P denoting the principal value. The sum of the imaginary part (identical to the modified 
zeta function) is rewritten as, 

_ gD/2 ^ ^D/2-1 

- ^"^^ = r^TiTTT^TTTTTTT M^:7-,^w^ • (3.3) 



n=l 



7r^/2-ir(D/2) Jo e*+W£ - 1 ' 
In D = 2, the integral can be performed explicitly, to give 

-^ln(l-e-'r) . (3.4) 

In view of (|3.3|) (or ( |3.4|) ), the imaginary part becomes significant when S > x [E > rn?), 
which is a rather well-known result from a one-body problem of the Dirac particle under 
an external potential — the Klein paradox. The appearance of the imaginary part implies 
the vacuum instability by means of external electric fields 0. The initial vacuum goes to a 
new one with emissions of particle pairs to neutralize the electric field which, however, has 
assumed constant so that the process will never end. To avoid the situation electric fields 
must be localized or the condition, 

8 <x , (3.5) 

must be assumed. (See Fig. |I|.) Moreover if a perturbative expansion in terms of the 
electric field is employed we cannot see any imaginary part at all since £^ = is an essential 
singularity. The program matches with our strategy. As was stated in the introduction, our 
interest lies in seeking a change when external fields are present, then it is enough to regard 
external fields to be very small: perturbations would be useful. 
Here we should make three comments. 



(i) According to the condition ( |3.5|) there is no need for an infrared cut-off in this case so 
that we put e = hereafter. 
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(ii) It is impossible to talk about criticality of electric fields , since it is defined through a 
massless condition under the variation of electric fields, which contradicts the condition 
T5D. 



(iii) The calculation is performed in the Euclidean world throughout. (In the above, we 
visited in the Minkowski space just because of glancing at the infrared structure.) 

With these spirits we first expand the integrand in such that 

tcotht = l + - + 0(t^) (3.6) 
3 

then integrate each term, and perform the analytic continuation £ —iS, to obtain 
D = 2 case: 



V2{x) = + xE,{-x) + - ^ 1 ^""^ - 1 ' (3.7) 



27r £^ 
— = f,{x) = Et{-x) + - , 



D = 3 case: 



2 ^i^)=E^i-x) + -i- + -,)e-^ , (3.8) 



V^{x) = -^x + 3 ( 4x - - j xV2r(l/2, x) + -(1 - 2x)e-- - - , (3.9) 



D = A case: 



4:71^ I x^ £'^\ 1 1 

^4(x) = - f Y - y ) E^hx) + -(1 - x)e-- - - , (3.11) 



- ^ = /4(x) = -xi?.(-x) - (l - Q . (3.12) 
Here use has been made of 

/CO roo 

dte-H-^ = dte-'T^ , (3.13) 

and 

dte-H-^'^ . (3.14) 

As was expected from the above argument there is no imaginary part. It should be noted 
that our perturbation expansion (|3.6|) is of course an asymptotic expansion to the effective 
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potential ( |2.26|) - (p.28|) since we have regarded t in ( p.6|) as t£ which apparently exceeds 



the radius of convergence when r — oo. 

We plot fnix) in Fig. |^ (a)-(c) for several fixed values of £. First note that in each plot, 
the thin-dashed line denoting fD{x)\£=x shows the lower bound of x (|3l5| ): the vacuum is 
unstable in the left region to the line. From the figures we can see that the external electric 
field makes the mass smaller for a given coupling. This can more easily be seen from Fig. |^ 
(a)-(c) where the relation between S and x are depicted for several fixed values of g"^. In 
each plot, the thin-dashed line denoting S = x, again shows the lower bound of x ( |3.5| ); the 
vacuum is unstable in the upper triangular region, g"^ in (b) and (c) are the critical couplings 
without any external fields in D = 3 and 4 and given by tt^/^/A and 47r^/A^ respectively. 

Finally it is necessary to discuss plausibility of our approximation. By changing the 
integration variable such as rx r , (|3.1|) is rewritten as 



X- dr T---^e-^ { — ) coth{ — ) . (3.15) 

Therefore our approximation becomes good when 

l^oo. (3.16) 

However, as is seen from Fig. ^ where the integral ( |3.1| ) for D = A with £ = 0.25 is depicted 
as the function of x, matching is excellent for any x down to £. 



IV. CONSTANT MAGNETIC FIELD 

The case for a constant magnetic field is considered in this section. (This reads, covari- 
antly in the Minkowski space, F^^F^^ > with F^^F^^ = in 4-dimension.) Contrary to 
the previous case, there arises no imaginary part in the effective potential: the integral in 



(|]2|) reads 



cirr---^e-"Vi3coth(ri3) (4.1) 

with B denoting a magnitude of the magnetic field in the A^ unit (0 < S < 1). There is 
no need for an analytic continuation so that any imaginary part does not occur. Therefore, 
unlike the electric case, there is no restriction on B and x (except for < S < 1 and 
< x < 1) such as ( |3.5| ), which enforce us to keep e non zero and to modify the naive 
expansion ( |3.6|) , since in this case B/x can become large contrary to the electric case where 
£ /x <1 from (p.5|). We thus arrange (|4.1|) such that 



dTT'-'^e-^'-'irB) coth(ri3) 

1 - e-2-« 



B r'drT-Te----+ /"drr-f-V^(^+2^)^^^ . (4.2) 



The last factor in the second integral can be expanded as 

t t 



e* - 1 2 12 



T^ + T7T + 0(n, (4.3) 
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and then all the term can be evaluated with the use of ( |3.13|) or ( |3.14|) . This gives an 
improved asymptotic expansion. The results are read as follows: 
D = 3 case: 



- . ^ 27r3/2 / e \ 1 



2{2x + B) + 



X + 2B 



+ -(1 -2x- S)e~(^-+26) _ 2B x^/^r{l/2, x) - e"^' 
3 



x + 2By/^r{l/2,x + 2B) 
2 



(4.4) 



27rV2 



2 - 



B 



B^ 



2 + 



X + 2B 6{x + 2By 
B^ 



[X 



+ 2BY/^r{l/2,x + 2B) 



3{x + 2B) 



-Bx-'/'r{l/2, 



X) 



D = 4 case: 



X 



A2 



-+xB + -\e,{-{x + 2B)) 



+ - x)e-("+2^) + B [xEi{-x) + e" 



1 
2 



(4.5) 



(4.6) 



h{x) 



[x + B)E,{~{x + 2B)) ^^fl^^^ e-(^+^^) + BE,{-x) . (4.7) 



We plot foix) in Fig. |^ (a) and (b) for several fixed values of B, from which we see that 
any point on the line fulfills the condition f']y{x*) > 0, ( p. 241) . The condition ( p.25| ) is 
also satisfied, since ( p.25| ) is always true for any monotonically increasing fnix). Since the 
minimum of x is e/A^, a finite critical coupling exists as a solution of ( p.23|) at x = e/A^; 



(47r) 



D/2 



(4.8) 



trlA^-2j^(,/A2^ 

for a fixed B. However, when e/A^ 0, /^(e/A^) behaves as 

Me/ A') ~ -S(e/A2)-V2r(i/2, e/A^) + O ((e/A^)) , (4.9) 
hie/ A') ~ BE,{-e/A') + O ((e/A^)) , (4.10) 

so that f£){e/A^) — oo, as far as S 7^ 0. This implies that the critical coupling goes to 
zero, Qc —* 0, for any non-zero B. The situation is similar to the 2-dimensional case without 
external fields, where from (|3l 



Mx) = E,i-x) 



(4.11) 



which is again divergent when x — 0. Keeping IR cut-off e finite, we obtain a finite critical 
coupling 
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which also becomes zero when t/h? 0. (Compare Fig. ^ (b) with the S = case of 
Fig. 1^ (a).) Gusynin, Miransky, and Shovkovy have interpreted this similarity in terms of 
the dimensional reduction 0. 

We can also see, from Fig. ^ (a) and (b), that the dynamically generated mass in a 
constant magnetic field is larger than that without any external field for a fixed coupling. 
This can more easily be seen from Fig. ^ (a) and (b), where the relations between B and x 
is depicted for several fixed values of g^. 

We adopt the improved expansion ([4.2|) with (|4.3| ) against the previous one (p.6|) . The 
expansion becomes exact when 1 +x/{2B) oo, which can be recognized by looking at the 
second term in the right-hand side of (|4.2| ). However as is seen from Fig. matching to the 
exact value is much more excellent than the previous case in Fig. ^, even if only the first 
three terms are taken into account. 



V. GENERAL CONSTANT FIELDS IN D=4 

So far, in D=4, the case F^^^F^'^ = (in the Minkowski metric), that is, B • E = 
has been assumed, but in this section a more general case F^^F^^ 7^ 0, is considered. We 
vary F^^F^'^ with F^^jF^'^ being fixed, which is interpreted, for example, such that the angle 
between B and E is shifted from 7r/2, while keeping the magnitude of B and E fixed. 

In this case, imaginary parts in the effective potential are also unavoidable. To see this, 
the analytic continuation \E\ — > — "il-El is performed in the third relation in ( p.l7|) and a 
covariant notation is adopted, such that 



F^^F^" = 2(B'-E 



F^^uF^" = -AB • E , (5.1) 



givmg 



(F^,FM-)2 + (F^,F/^-)2 + F^^F^^^ = F+ , (5.2) 



F^ -sgn(F^,F'^^) V V(^M-^'^'')' + {F^uFf^-^y - F^,F>^'^ = -iF_ , (5.3) 



where sgn(x) designates the sign function 

^^■'(^■) - { -\ t < 0) : <^-'' 

and the notation F^ has been employed to distinguish the Euclidean quantities from the 
Minkowski ones. Since there exists a suitable Lorentz frame where \B'\ = \E'\ = \F-\, 
and B' \\ E' ,we regard F^ {F-) as a Lorentz invariant magnetic (electric) field. The integral 
in (|2.28|) becomes, after the analytic continuation, to 
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/oo 



P / rfrr"''e"™'(^^-^+-^-)coth(rJ^+)cot(rJ^_) 



-tTT > 



(5.5) 



n=l 

with = and = \F^\/A'^, obeying < J^± < 1. (The sign of F„ or F^^F'^'' 

is irrelevant because of parity invariance.) From (|5.5| ), JF_ is solely responsible for the 
imaginary part as is expected. (When J-"^ there remains the imaginary part.) Therefore, 
as was done in the electric field case, we must set the condition 



(5.6) 



for not having a large imaginary part, which enables us to put e — > and again prevent us 
from talking about criticality of the coupling and external fields. 

Now go back to the Euclidean world and regard and J^^ as perturbations as before. 
Follow the same procedure as in the preceding section: first expand the integrand in the 
left-hand side of (|5.5|) such that 



/oo 

"(rj^^) coth(rJ^^) 



dTT ^e ™ 



'1 

/oo 
dr r" 



•3g-r(z+2.Ff) 



r^^ + ^(r^+ 



E\2 



{tJ^^) coth(rJ^^ 



(5.7) 



where we have employed the relations ( [4. 21) and (|4.3|) developed in the pure magnetic field 
case, which is an extremely good expansion. Then we expand {tJ-'^) coth(rjF:^) by using 
the expansion (|3.6| ), integrate each term, and make analytic continuation, J-'^ —ij-'^, to 
obtain 



"1 

+ 



xEii—x) +1 e 



X 



1 + 



X _ 
1 \ 



3 



2 3(x + 2J'+) 9(x + 2J^+) ' x + 2J^_ 
As was expected, there is no imaginary part. The gap equation is then found as 
47r2 



Ei{-{x\2T^)) 
\ ■ (5.8) 



= Uix) 

T'^ ( 1 
Fz(-x) + ^ l + -)e^ 
3a; V X, 



^ 3(x + 2J^4 



(x + J^+)Fi(-(x + 2J^+)) 
1 \ 



9(x + 



1 



3(a; + 2J^+) 
2 



X + 2J^4 



2J^4 



2^,) 2 



(5.9) 
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The results for the cases, Tf.yT'"' = 0.2 (magnetic-hke) , J^^^J^^'^ = 0, and J^^^J^'^^ = -0.2 
(electric-hke) are shown in Fig. ^ (a)-(c), from which we first notice that the dynamical mass 
becomes always smaller as the magnitude of jF^j^jF^^^ goes larger. This could be understood 
from the following facts: (i) from ( ^.3|) , JF_ is a monotonically increasing function of \ J^^u^'^''\ 
for a fixed value of JF^i^JF^*^. (ii) JF_ is nothing but covariant electric field so that JF„ becomes 
larger the mass goes smaller according to the discussion on the pure electric field case in Sec. 
[in| . In this respect it should be noted that there is almost no role of J-'+, also an increasing 
function of \J-'fj,u^'^'^\, in this phenomena. This might be seen from the form of the imaginary 
part of the effective potential (p.5[), where J-'_ is essential. Second, we notice that in Fig. 8 
(a) ((c)) graphs in the physical region defined by the condition ( |5.6| ) are shifted to the larger 
(smaller) mass side with respect to the curve of no external field, jF^jyjF^*^ = J^^^^^" = 0, 
which reflects the fact that in the magnetic-like case, T^u^^^ > 0, the mass goes larger 
while in the electric-like case, T^yJ^^^ < 0, it goes smaller. The thin-dashed line indicates 
the lower bound of x ( |5.6| ). 

It should also be noted that the asymptotic expansion in this case matches excellently 
with the exact value for almost all values of x down to JF_. (See Fig. P.) 



VI. CONCLUSION 

We examine the NJL model in 2-, 3-, and 4-dimension coupled to external constant 
electromagnetic fields. The results are summarized as follows: an electric (magnetic) field 
reduces (raises) the dynamically generated masses, that is, electric fields oppose xSB, while 
magnetic fields enhance it. 

In the case for a pure magnetic field (in 3- and 4-dimension), we obtain a well-defined 
effective potential with a UV cut-off as well as an IR cut-off e. When e is nonzero, we 
obtain a non-zero critical coupling which, however, goes to zero as e ^ 0. We thus reconfirm 
that the critical coupling is zero for any infinitesimal magnetic field. 

We find that the effective potential has the imaginary part in the electric case defined by 
condition F^^F'^" < in 2- and 3-dimension and F^^F^^ < or F^^F^^'^ 7^ in 4-dimension. 
In these cases, imaginary parts should be small and negligible, so that the condition, £ < x 
as well as JF_ < a; in 4-dimension, has been obtained, which enables us to put the IR cut-off 
e zero. Moreover, this condition prevents us from selecting out the critical quantities, such 
as critical couplings or critical fields, since those are defined through the transition from a 
massive to a massless state or vice versa. 

However the imaginary part cannot be seen under the perturbation theory with respect 
to the external fields. This perturbative expansion, realized as an asymptotic expansion, 
excellently matches with the exact value by adopting only the first few terms. We have 
assumed that the external fields are constant, so that we can calculate the functional deter- 
minant exactly. If, however, the fields depend on space-time, we cannot calculate it without 
approximation. In a usual approach, a weak field approximation is employed, which would 
therefore be a fairly good expansion according to our analysis. 

In our present work, we can ignore the chiral anomaly, which is trivial in the abelian 
case but indispensable to investigate the reality of the dynamics of QCD in the low energy 
region. It is also necessary to work with non-constant gauge field, such as the instanton 
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configuration. Generalization of tlie present work to a non-abelian gauge field and including 
the effect of the anomaly is our next step. 



APPENDIX A: CALCULATION OF Gd{tF) 

In this appendix the calculation of Gd{F) with a real antisymmetric tensor F^^ is ex- 
plicitly given. 

The definition for Gjj^F) is 



Gd{F) 



tr(e^°"'"'-'^f"'^ 
trl 



det 



sin F 



(Al) 



where tr is for gamma matrices and det is for F^v a.s a D x D matrix. The result is 



Gd{F) 



Focoth Fd 



for D = 2,3 



F+F_ cothF+ cothF_ for = 4 , 



(A2) 



with 



F2 = E , F^ = ^B^ + E' , 
F± = {\B + E\±\B - E\} /2 



(A3) 



where E (or E) and B (or B) denote an electric field and a magnetic field respectively. 

The derivations of (|A2|) and ( |A3|) are shown as follows: m D = 2 case, they can easily 
be obtained with 



7m = 



M 1 



75 



(A4) 



and with 



E 
-E 



(A5) 



which can be diagonalized into i-E'diag[l, —1]. Therefore 



G2{F) 



tr e 



trl 



det 



/ sinhi?' 



EcothE 



In = 3 case, since F^u in the 3-dimensional space-time is expressed as 



(A6) 



/ B El 

-B ^2 I = ii^E^Ti - E1T2 + BTs) 
\-Ei -E2 



(AT) 



where Tj(i = 1, 2, 3) is a basis of the SO (3) algebra in the adjoint representation, it can be 



diagonalized into iy B'^ + E diag[l, 0,-1]. Therefore 
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The calculation of the trace part of ( [A^ ) can be done with the use of gamma matrices 
in 3- dimension, given by 



7m 



Therefore 



_^ j, (Tfiu = ^[l^„lu] = ef,up(Tp; /i, i/,p = 1,2,3. 







(A9) 



1 F -If ( \ _ E2cri - Eia2 + Ba-s ^ , , . 



which can be diagonalized as y B'^ + £^^diag[l, —1, 1, —1], yielding to 

\iT( G^^^^'^ til^^ 



trl 



cosh F'- 



3 • 



Combining ([AED and ( |A11|) , we obtain and (|X3D in D = 3 case. 
In D = 4 case, F^,^ is expressed as 



( 53 -52 El ^^ 

— i?3 i?i ^2 
i?2 — Bi 

V -El -^2 -^3 y 



i{B ■ M + E- N) 



pv 



(All) 



(A12) 



where Mj and iVj (i=l,2,3) is a basis of the 5*0(4) algebra in the fundamental representation. 
Since Mj and Ni are decomposed into the subalgebra 



satisfying 



[Jj, Jj] = ieijkJk , [Ki, Kj] = ieijkKk , [Jj, -ft'j] = 



(A13) 



(A14) 



F^y, now expressed as i {{B + E) • J + {B — E) • K}^^, can be diagonalized into 
idiag[F+, — F_, —F_]. Therefore 



det 



/ sinF 



F.F^ 



sinh Fj. sinh F_ 



(A15) 



As for the trace, in order to diagonalize anyFav, we first define 



M^, = ^e^upCTup , J^^L = ]^(ypA ; /i, z/, p = 1, 2, 3 . 



(A16) 



Similar to M and A/", M.i and Mi can be decomposed into the subalgebra. 
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which also satisfies the same algebra as ( |A14|) . Thus 



(A17) 



(A18) 



which can be diagonalized into diag[|B + E\^ —\B + \B — E\^ —\B + E\]. Hence 



trfea'^'"'^'"') 1 , , 1 , 

— = - cosh \B + E\ + - cosh \B — E\ = coshF. coshF_ 

trl 2 ' ' 2 ' ' + 



(A19) 



Combining (|A15D and (|A19|) , we finally obtain (|A2D and ( |A3|) for D = 4 case 
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FIGURES 

0.1 0.2 0.3 0.4 0.5 




FIG. 1. Behavior of the imaginary part of the potential: we write the quantity of (|3.3| ) as 
—ImVDix) and put £ = 0.1. The imaginary part is exponentially dumping as x goes larger and is 
smaller than 10~^ when x > £. ImVDiO) is finite in D = 3 and D = A but divergent in D = 2. 




FIG. 2. foix) for a pure electric field in (a) D = 2^ ih) D = 3, and (c) D = 4 for several fixed 
values of £. In each plot, the thin-dashed line denotes fD{x)\£=x and shows the lower bound of x 
(|3.5|): the vacuum is unstable in the left region to the line. 
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0.2 ■ 0.2 




FIG. 3. Relations between the magnitude of the electric field £ and the dynamically generated 
mass X for several fixed values of in (a) D = 2, (b) D = 3, and (c) D = 4. In each plot, the 
thin-dashed line shows the lower bound of x ( |3.5| ). in (b) and (c) denote the critical couplings 
without external fields in each dimension and are given by vr^/^/A and 47r^/A^ respectively. This 
shows that the mass becomes smaller when the electric field goes larger. 

Ir , -0.01 



0.8 ; - 0.008 
x/£»l 




I . . , , , 

0.2 0.4 0.6 0.8 1 
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FIG. 4. Plausibility of the asymptotic expansion in the pure electric field case: we write the 
integral (|3.lD in = 4 as J{x) = £ J^dr T~'^e~'^^ coth(rf ) and put £ = 0.25. The thin line denotes 
the (numerically evaluated) exact value and the dashed line represents the asymptotic expansion up 
to the second term. The dotted line denotes the relative difference [{Approx.) — (Exact)]/ {Exact) 
using right-hand scale. Matching is excellent for any x down to £. 





(b) 

FIG. 5. foix) for a pure magnetic field in (a) D = 3 and (b) D = A for several fixed values of 
B. The minimum of x is e/ so that there exists a finite critical coupling as far as e 7^ 0. However 
when e ^ critical couplings go to zero (foix) —00) for any non-zero B. 



si'=Q.SgJ' sf=g,' 9^=1.6 9^=0,7 9;^ 9==9,^ 9^=1-3 9o' 




FIG. 6. Relations between the magnitude of the magnetic field B and the dynamically generated 
mass X for several fixed values of in (a) -D = 3 and (b) D = 4. g'^ denotes the critical couplings 
without external fields given by vr^/^/A in 3-dimension and 47r^/A^ in 4-dimension. This shows 
that a magnetic field prompts the generation of mass thus enhances xSB. 
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0.2 0.1 0.6 0.8 1 



FIG. 7. Plausibility of the improved asymptotic expansion in the pure magnetic field case: we 
write the integral in D = 4 as J(x) = B /^dr r-^e"™ coth(TS) and put B = 0.5. The thin 
line denotes the (numerically evaluated) exact value and the dashed line represents the improved 
asymptotic expansion up to the first three terms. The dotted line denotes the relative difference 
[{Approx.) — [Exact)]/ {Exact). Matching is significantly excellent for all values of x. 



m fix) 




X 
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FIG. 



0. 



/4(^) 

and (c 



for several fixed values of T^y!F^^^ for the cases, (a) T^j,i,T^^ = 0.2, (b) 



-0.2. The thick hne, T^yT^''' 



0, that is, no ex- 



ternal fields, is also depicted in comparison. The thin dashed line indicates the lower bound of 
X (^^), which is obtained by putting = x and Tj^ = \Jx'^ + ^^u^^^ /'^ in fiix). These show 

that J-p^uJ-^'^ opposes mass generation. In (a) ((c)), graphs in the physical region, satisfied with 
( |5.6| ), are shifted to right (left) comparing to the curves in (b). This reflects the fact that in the 
magnetic-like case, T^yT^^ > 0, the mass goes larger while in the electric like case, ^F^j.^J^'^'^ < 0, 
it goes smaller. 




CD 
< 



FIG. 9. Plausibility of the asymptotic expansion for general external fields in D 



4: we write 
coth(rjF_) and put JF^ = 0.4 



the left-hand side of (5.5) as J{x) = T^T- /^°°dT e~™ coth(rjF 
and = 0.2. The thin line denotes the (numerically evaluated) exact value and the dashed line 
represents the asymptotic expansion in the same approximation as the one for (^. j). The dotted 
line denotes the relative difference [{Approx.) — {Exact)]/ [Exact). Matching is excellent for almost 
all values of x down to !F- . 
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